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Abstract 

The frequency of inherited bilateral autosomal recessive non-syndromic hearing loss (ARNSHL) in Pakistan is 1.6/1000 
individuals. More than 50% of the families carry mutations in GJB2 while mutations in MY015A account for about 5% of 
recessive deafness. In the present study a cohort of 30 ARNSHL families was initially screened for mutations in GJB2 and 
MY015A. Homozygosity mapping was performed by employing whole genome single nucleotide polymorphism (SNP) 
genotyping in the families that did not carry mutations in GJB2 or MY015A. IVlutation analysis was performed for the known 
ARNSHL genes present in the homozygous regions to determine the causative mutations. This allowed the identification of 
a causative mutation in all the 30 families including 9 novel mutations, which were identified in 9 different families {GJB2 
(c.598G>A, p.Gly200Arg); MY015A (c.9948G>A, p.Gln3316Gln; c.3866+lG>A; C.87670T, p.Arg2923* and c.8222T>C, 
p.Phe2741Ser), TMCl (c.362+18A>G), BSND {c.97G>C, p.Val33Leu), TMPRSS3 (C.7260G, p.Cys242Trp) and MSRB3 (c.20T>G, 
p.Leu7Arg)). Furthermore, 12 recurrent mutations were detected in 21 other families. The 21 identified mutations included 
10 (48%) missense changes, 4 (19%) nonsense mutations, 3 (14%) intronic mutations, 2 (9%) splice site mutations and 2 (9%) 
frameshift mutations. GJB2 accounted for 53% of the families, while mutations in MY015A were the second most frequent 
(13%) cause of ARNSHL in these 30 families. The identification of novel as well as recurrent mutations in the present study 
increases the spectrum of mutations in known deafness genes which could lead to the identification of novel founder 
mutations and population specific mutated deafness genes causative of ARNSHL. These results provide detailed genetic 
information that has potential diagnostic implication in the establishment of cost-efficient allele-specific analysis of 
frequently occurring variants in combination with other reported mutations in Pakistani populations. 
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Introduction 

Deafness or hearing loss is tlie most common congenital 
sensorineural disorder afTecting about 1 in 1000 children. Genetic 
factors contribute to approximately half of the cases of hearing loss 
[1], which can be either syndromic or non-syndromic. The former 
is responsible for about 30% of prelingual deafness in combination 
with abnormalities of other organs. Non-syndromic deafness is 
usually due to abnormalities of the middle and/ or the inner ear 
and is found in 70% of the hereditary cases [2]. The disease is 
genetically heterogeneous and currendy 134 loci have been 



determined and 80 genes identified for the non-syndromic type 
(Hereditary Hearing Loss Homepage, URL: http: / / 
hereditaryhearingloss.org/). More than 400 syndromes are known 
with deafness as one of the symptoms and for many of these the 
causative genes have been identified (Online Mendelian Inheri- 
tance in Man: http://www.ncbi.nlm.nih.gov/omim). 

Syndromic and non-syndromic hearing loss display an autoso- 
mal recessive, autosomal dominant, X-linked, Y-linked or 
mitochondrial pattern of inheritance [3]. Eighty percent of the 
early onset non-syndromic cases have an autosomal recessive 
inheritance pattern (ARNSHL), while autosomal dominant. 
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X-linked and mitochondrial non-syndromic deafness contribute to 
18%, 1-3% and <1% of the cases, respectively [4]. In Pakistan, 
hearing impairment is severe and congenital in 70% of the cases 
and the increased occurrence of these conditions is due to a high 
rate of consanguineous marriages (60%); profound bilateral 
deafness occurs at 1.6 per 1000 iiidi\ iduals [5]. 

The genes most frequendy involved in ARNSHL are those 
encoding gap junction protein beta 2 {GJB2, MIM# 121011), 
myosin XVA (MYOISA, MIM# 602666), transmembrane chan- 
nel-like 1 {TMCl, MIM# 606706), solute carrier family 26 (anion 
exchanger) member 4 (SLC26A4, MIM# 605646), otofi-rlin 
{OTOF, MIM# 603681) and cadherin-related 23 {CDH23, MIM# 
605516), each of which has been found to contain more than 20 
different mutations, most of which have been reported in 
consanguineous families [6]. Mutations in GJB2 are the most 
common cause of ARNSHL and explain up to 50'M) cases in the 
Mediterranean regions [7,8,9] while mutations in multiple genes 
have been shown to cause deafness in the remaining cases. In the 
Pakistani population mutations in MT015A account for 5% of the 
recessive deafness [10]. 

In the present study a panel of 30 unrelated consanguineous 
Pakistani families was initially tested for involvement of GJB2 and 
MY015A followed by whole genome homozygosity mapping and 
candidate gene sequencing. This approach resulted in defining the 
mutation spectrum of the disease in the current panel. 

Materials and Methods 

Ethics statement 

The current study conformed to the tenets of the Helsinki 
declaration and was approved by the Department of Biosciences 
Ethics Review Board of the COMSATS Institute of Information 
Technology, Islamabad, Pakistan. All patients, their normal 
hearing family members and 89 ethnically matched control 
individuals were informed about tlu' purpos(; of the study and 
written consent was taken before recruitment and sampling. 
Informed written consent of minors was obtained from their 
guardians. 

Genotyping 

A total of 30 consanguineous families with ARNSHL were 
ascertained from different regions of Punjab, Pakistan. Audiom- 
etery was performed on a few members from each family to 
determine the level of hearing loss. 

Blood samples were collected in EDTA containing vacutainers. 
DNA extraction from these samples was carried out using a 
standard phenol-chloroform/organic method [11]. Microsatellite 
markers were analyzed as described previously [12] and Sanger 
sequencing was performed according to Schraders et al. [13]. 
Primers were designed using Primer3 software (URL: http://www. 
bioinformatics.nl/cgi-bin/primer3plus/primer3plus.cgi/; [14]). Prim- 
er sequences are a\ ailable upon request. PGR was performed using 
standard conditions. 

Prior to whole genome single nucleotide polymorphism (SNP) 
mapping all families were pre-screened for mutations in (JJB2 by 
Sanger sequence analysis. For the MYOISA locus, microsatellite 
markers (D17S1843, D17S2196, D17S783 and D17S1824) were 
genotyped and in families in which haplotype analysis showed 
compatibility with genetic linkage, Sanger sequence analysis was 
performed for MYOISA. For family 1 IDF, microsatellite markers 
(D9S166, D9S1806 and D9S1876) were used for exclusion of the 
TMCl region. 

Families with no identified mutation in (J}B2 and MYOISA 
were further genotyped using the lUumina HumanOmniExpress 



whole genome single nucleotide polymorphism (SNP) microarray 
(>700 K SNPs) or the IHumina Human Linkage- 12 panel 
according to the manufacturer's protocols. Homozygosity map- 
ping was performed using an online tool Homozygosity Mapper 
(URL: http://www.homozygositymapper.org; [15]). The loga- 
rithm of odds (LOD) score for the homozygous regions identified 
in family DFR18 was calculated using the Gene Hunter v2.1r5 
program in the easyLINKAGE plus v5.08 software package [16]. 
The exons and exon-intron boundaries of the known candidate 
genes {BSND (NM_057 176.2), GJB2 (NJvI_004004.5), HGF 
(N]\'I_000601), MYOISA (NM_()16239.3), MSRB3 (NM_00 103 1679.2), 
SLC26A4 (NM_000441.1), TMCl (NM_138691.2), TMPRSS3 
(NM_024O22.2) and TMIE (NM_147196)) present in the homozy- 
gous regions were sequenced in the proband of the families. 
Segregation analyses for identified mutations in the correspond- 
ing families were performed by Sanger sequencing except for 
family DFR24 and DFR18. In these families, the segregation of 
TMPRSS3 and MSRB3 mutations was analyzed with restriction 
digestion with the enzymes Acil and Tsel, respectively (New 
England Biolabs Inc. UK). 

Eighty nine ethnically matched controls were sequenced for the 
novel mutations identified in the current study. Control panel 
screening of GJB2 (c.598G>A); MYOISA (c.9948G>A), (c.3866+ 
1G>A), (c.8767G>T), (c.8222 T>C); TMCl (c.362+18A>G) and 
BSND (c.97G>G), variants was done by Sanger sequencing. While 
for TMPRSS3 (c.726C>G) and MSRB3 (c.20T>G) variants Acil 
and Tsel PCR-RFLP were performed, respectively. 

In silico prediction of the identified variants 

In silico prediction of the identified variants was performed using 
online prediction tools. Sorting intolerant from tolerant (SIFT: 
http://sift.jcvi.org; [17]) and Polymorphism Phenotyping v2 
(PolyPhen-2: http://genetics.bwh.harvard.edu/pph2; [18]) were 
used for analyses of missense changes. Effects on splicing were 
evaluated with NetGene2 (http://www.cbs.dtu.dk/services/ 
NetGene2/; [19]). SignalP 4.0 (http://www.cbs.dtu.dk/services/ 
SignalP-4.0/) was employed to predict the presence and location 
of the signal peptide /non-signal peptide cleavage sites [20]. In 
addition, the online tool. Have your Protein Explained (HOPE) 
(http://www.cmbi.ru.nl/hope/input) was used to predict the three 
dimensional structural changes at the protein level [21]. The 
exome variant server (EVS) and an in-house exome database 
(Human Genetics, Radboud University Medical Gentre) were also 
searched for the presence of putative pathogenic variants. 

Minigene Construction and Splicing Assay (for C.362+ 
18A>G mutation in TMCl and c.9948G>A mutation in 
MYOISA) 

A plasmid containing the genomic region encompassing exons 
3-5 of RHO inserted at the EcoRI/Sall sites in the pCI-NEO 
vector was used for in vivo splicing assays [22]. The plasmid was 
adapted from the Gateway cloning technology (Life technologies) 
according to the manufacturer's protocol. PGR amplified 
fragments of wild-type and mutant TMCl exon 8, along with 
flanking intronic sequences were generated with the following 
primers, 5'-GGGGAGAAGTTTGTAGAAAAAAGCAGGGTT- 
Gtgggtcctaatgttgactgc-3' and 5'-GGGGAGGAGTTTGTAGAA- 
GAAAGCTGGGTGggatttagaaaatcaatatcaggg-3', that contain 
the attBl and attB2 sites necessary for Gateway cloning. Similarly, 
amplified fragments of wild-type and mutant MYOISA exon 61 
were generated using primers, 5'- GGGGAGAAGTTTGTA- 
GAAAAAAGGAGGGTTGttccccaggagaaatggag-3' and 5'- GGG- 
GAGGAGTTTGTAGAAGAAAGGTGGGTGcagcttgctgggaag- 
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Table 1. Spectrum of recurrent GJB2 mutations in Pakistani families with autosomal recessive non-syndromic hearing loss 
(ARNSHL). 





Mutation identified (Protein change) 


Type of mutation 


No. of 
families 


No. of affected 
members 


Frequency in EVS 


C.710A (p.Trp24») 


Nonsense (homozygous) 


5 


26 


Absent 


C.2310A (p.Trp77*) 


Nonsense (homozygous) 


4 


14 


AA = 0/AG = 1/GG = 4299 


c.35delG {p.Glyl2Valfs»2) 


Frameshift (homozygous) 


2 


8 


Absent 


c.35delG (p.Glyl2Valfs*2) c.439G>A (p.Glul47Lys) 


Missense (compound heterozygous) 


1 


3 


Absent 


c.380G>A (p.Arg127His)* 


Missense (homozygous) 


2 


4 


AA = 0/AG = 26/GG = 4274 


c.377_378insATGCGGA (p.Arg 1 27Cysfs*85) 


Frameshift (homozygous) 


1 


2 


Absent 



As reference sequence NM_004004.5 was employed. EVS, exome variant server; 
*The pathogenicity of this mutation is controversial. 
doi:10.1371/journal.pone.0100146.t001 



gac-3'. Recombinant vectors were employed for transfection of 
HEK293T cells that were incubated at 37°C for 24 hours. RNA 
was then isolated from the transfected cells and reverse 
transcribed into cDNA and sequenced. Forward primer 5'- 
cggaggtcaacaacgagtct-3' and reverse primer 5'-aggtgtaggggatgg- 
gagac-3' were used to amplify and sequence the amplified 
cDNA fragments along with flanking RHO sequences as present 
in the vector. 

Results 

AH families in this study were consanguineous (Table SI) and all 
patients in these families were diagnosed with severe to profound 
congenital hearing loss. Causal mutations were identified in all 30 
families and the 9 novel mutations identified in the current study 
were not present in any of the 89 ethnically matched control 
individuals. 

Sequence analysis of GJB2 in the current cohort identified 16 
families with mutations in this gene (Table 1 ; Table 2) segregating 
with hearing loss. A recurrent nonsense mutation c.71G>A 
(p.Trp24*) was the most common and was found in 5 families 
(31%), followed by another nonsense variant, c.231G>A 
(p.Trp77*) that was identified in 4 families (25%; Table 1). In 
addition to the 6 recurrent GJB2 mutations, a novel homozygous 
change c.598G>A (p.Gly200Arg) was found in one consanguin- 
eous family DFRIO (Table 2); as predicted by the HOPE server 
the larger side chain of the mutant residue arginine might well 
affect the proper folding of the cysteine rich domain. In addition, 
the residue resides in the conserved region of the protein (Figure 1 ; 
Figure SI A). 

Out of the 9 novel mutations (Table 2) identified in the current 
study, 4 were present in MT015A, of which 2 were splice site 
mutations (c.9948G>A and c.3866+lG>A), one was a nonsense 
(c.8767C>T) and one a missense mutation (c.8222 T>C). The 
c.9948G>A variant changes the last nucleotide of exon 61 and is 
predicted to affect the splice donor site. The splice donor site in the 
reference sequence has a highly confident score of 94%, which is 
reduced to a score of 24% in the mutant. To determine the effect 
of the c.9948G>A mutation on splicing, a minigene approach was 
used. This showed correct sphcing of the wildtype MT015A exon 
61, while the c.9948G>A mutation almost completely abolished 
the normal sphcing (Figure 2). NetGene2 predicted the abolition of 
the splice donor site of exon 5 as a result of the c.3866+lG>A 
mutation in intron 5 of MT015A. The nonsense mutation 
c.8767C>T (p.Arg2923*) was novel and predicted to lead to the 
synthesis of a truncated protein, while the missense mutation 



c.8222 T>C, leads to the substitution of serine for phenylalanine 
(p.Phe2741Ser) at amino acid position 2741 that resides in the 
conserved region of the protein (Figure 1; Figure SIB). This 
missense change was predicted to be deleterious by SIFT and 
Polyphen2 (Table 2). 

Haplotype analysis of STR-markers flanking TMCl showed 
compatibility with genetic linkage for family 1 IDF. Since sequence 
analysis of TMCl revealed a novel intronic mutation, c.362+18A> 
G, this family was not further analyzed by SNP-array genotyping. 
The variant segregated with the hearing loss in the family (Figure 1) 
and was predicted to create a novel splice donor site with a similar 
confidence score as the original splice donor site. To determine the 
effect of the c.362+18A>G mutation on splicing, a minigene 
approach was used. This revealed correct splicing of the wild-type 
TMCl exon 8, while the c.362+18A>G mutation resulted in a 
1 7 bp extension of exon 8 (Figure 3). This leads to a frameshift and 
a premature stop codon (p.Glul22Tyrfs*10). 

A novel homozygous missense mutation was identified in 
TMPRSS3 (MIM_605511), C.7260G (p.Cys242Trp) that co- 
segregated with hearing impairment in family DFR24 (Figure 1) 
and was predicted to be probably damaging by Polyphen2 and 
deleterious by SIFT (Table 2). This mutation is located in the 
peptidase S 1 domain and the HOPE server predicted the abolition 
of the catalytic activity of TMPRSS3. In addition the amino acid 
substitution is present in a region of the protein that is conserved 
across different species (Figure S 1 D) and therefore probably affects 
the core structure of the peptidase domain (Figure 4). 

In family DFR18 SNP microarray data analysis revealed a 
17.6 Mb homozygous region on chromosome 12 flanked by SNPs 
rs7978381 and rs7976686, a LOD score of 2.7 was calculated for 
this region. The known deafness gene, MSRB3 (MIM_613719), in 
the region was subsequendy sequenced, revealing a novel 
homozygous nucleotide substitution c.20T>G (p.Leu7Arg) in 
exon 4 (Figure 1; Table 2). The leucine at position 7 is located in 
the signal peptide of the MSRB3 protein and as a result of the 
substitution by arginine this signal peptide loses its function as 
predicted by SignalP 4.0. Furthermore, the Leu7 residue is 
conserved across species (Figure SIE). The mutation abolishes a 
7^f«/ restriction site, which allowed the segregation of this variant in 
the family to be checked by restriction digestion. By using the same 
analysis, this variant was also found heterozygously in 3 out of 1 78 
ethnically matched control alleles. 

The novel BSMD (MIM_606412) missense mutation c.97G>C 
(p.Val33Leu) in family 7DF (Figure 1) was predicted to be 
deleterious by SIFT and Polyphen2 (Table 2) and Val33 is 
conserved across species (Figure SIC). 
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Figure 1. Pedigrees and tiie segregation of novel mutations in l<nown deafness genes. Unfilled circles indicate unaffected females, unfilled 
squares indicate unaffected males, filled circles indicate affected females, filled squares indicate affected males, double lines represent 
consanguineous marriages, slashed line across the symbols indicate deceased individual, + indicates wild type allele, M indicates mutant allele. 
doi:l 0.1 371 /journal.pone.01 00146.g001 
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M = 100 bp marker 

1 =wildtype 

2 = mutant 
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Figure 2. Effect of IV1Y015A c. 9948G>A using a minigene approach. An agarose gel containing RT-PCR products detected from HEK293T 
cells transfected with the wildtype and mutant minigene construct and a schematic representation of the identified splicing products. The RT-PCR 
products were verified by sequence analysis. The C.99480A mutation leads to skipping of exon 61. 
doi:10.1371/journal.pone.0100146.g002 

Recurrent mutations in TMC; (c. 11 14G> A; p. Val372Met [23] (c.l337A>G; p.Gln446Arg [26]) and TMIE (c.241C>T; 
and c.lOOOT; p.Arg34* [24]), //GF(c.482-H991_2000delGAT- p.ArgSlCys [27]) also segregated in the respective families and 
GATGAAA and c.482-H986_1988delTGA [25]), SLC26A4 most likely are the disease causing mutations (Table 3). 



AT GG AGAAG AAAATT GAAGTTCTCAAT T CCGG AACTGCATGC TCAC CAC C 
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TTCCGGAACT 














< 




> 





n ryWCf/lntronSi' 
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Figure 3. Effect of TMC1 intronic mutation c.362+18A>G using a minigene approach. Electropherogram of the partial cDNA sequence of 
RNA derived from cells transfected with the pCI-NEO with either the mutant or wildtype JMC1 exon 8. The mutation leads to the insertion of 1 7bp at 
the 3'end of exon 8, which can be predicted to result in a premature stop codon in exon 9 (p.Glu122Tyrfs*10). 
doi:1 0.1 371/journal.pone.01 001 46.g003 
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Figure 4. Predicted effect of mutation C.7260G (p.Cys242Trp) on the tKiree dimensional structure of TiVIPRSSS. A) Wild type protein 
structure with an intact disulphide bridge showing position of the mutated residue (magenta). B) Close-up view of the structure showing the wild 
type residue cysteine (green) and the mutant residue tryptophan (red). In case of the mutant residue there will be no disulphide bridge at this 
position. 

doi:10.1371/journal.pone.0100146.g004 



Discussion 

In a cohort of 30 ARNSHL families the genetic defects were 
identified in 9 known deafness genes; GJB2, MTOISA, TMCl, 
BSND, TMPRSS3, MSRB3, HGF, SLC26A4 and TMIE. In the 
current panel recurrent as well as novel mutations were detected, 
the novel mutations were identified in GJB2, MY015A, TMCl, 
BSND, TMPRSS3 and MSRB3. 

GJB2 was the most frequently mutated gene in these families. In 
this gene the novel mutation c.598G>A (p.Gly200Arg) affects a 
residue in the cysteine-rich domain that is involved in the 
formation of intramolecular disulphide bonds [28,29]. Although 
the effect of glycine on the intramolecular disulphide bond 
formation cannot be predicted, the mutant residue arginine may 
affect the proper folding of the cysteine-rich domain because of the 
larger size of the arginine side chain. 

In two famihes the GJB2 variant c.380G>A (p.Argl27His) was 
found to segregate with the hearing loss, being present homozy- 
gously in 4 affected members. Other studies have also reported 
deafness patients carrying this mutation homozygously [29,30,3 1] . 
Based upon the higher carrier frequency in patients as well as 
controls, Padma et al. [30] suggested that it is unlikely that this 
variant plays a pathogenic role in deafness. However, Matos et al. 
[32] have proposed that this mutation is likely to be pathogenic 
when influenced by other unknown genetic or environmental 
factors [32]. In the EVS, however, no homozygous occurrence of 
the mutated allele has been reported to date (Table SI). 
Identification of 4 p.Argl27His homozygous patients out of 125 
deafness patients in the current study and the previously reported 
cases of 9 heterozygous cases and 1 compound heterozygous case 
in 70 deaf patients by Bukhari et al. [31], demonstrates that this is 
the most frequently occurring mutation in Pakistani deafness 
families and therefore becomes important in genetic counseling of 
Pakistani deaf patients. 

In the current study GJB2 mutations were found to be the most 
common, followed by MTOISA mutations. The overall frequency 
of the two most common nonsense variants c.71G>A (p.Trp24*) 
and c.231G>A (p.Trp77*) of GJB2 in Pakistani deafness patients 
was obtained from the current (n=125) as well as previously 
reported studies of Santos et al. [33] (n = 430) and Bukhari et al. 
[3 1] (n = 70). Using these data, the frequency was found to be 
5.9% (p.Trp24*) and 4.3% (p.Trp77*) [31,33]. The variant 
p.Trp24* has a high prevalence in the Indian population as well 
[34], Pakistan and India have a shared genetic ancestry [35], 



which could be the reason of the presence of this variant in both 
populations. The third mutation p.Glyl2Valfs*2 identified in the 
current study has previously been reported from Northern areas of 
Pakistan [31] and also in Caucasians and Turks [36], thus 
indicating a possible founder effect of this mutation. In the current 
cohort, 35% (60/173) individuals were found to be carriers of 
GJB2 mutations. Therefore, as a preliminary screening of deaf 
families the sequencing of this gene is suggested in Pakistan. In the 
current study MTOISA was found to be the second leading cause 
of deafness in the Pakistani population, but still no recurrent 
mutation was identified in this gene. 

A total of 7 nonsense mutations were identified in GJB2, 
MTOISA and TMCl, these nonsense mutations are likely to cause 
nonsense mediated decay (NMD) because they are either present 
in the middle or near the 5 'end of the gene. 

In MTOISA two splice site, a nonsense and a missense mutations 
were identified. The c.9948G>A (p.Gln3316Gln) mutation affects 
the last nucleotide of exon 6 1 and changes the consensus splice site 
sequence. Based on the results of the minigene assay the 
c.9948G>A mutation is expected to lead to skipping of exon 61, 
which would result in a frameshift and NMD can be expected to 
occur at least for part of the mRNAs. The second splice site 
mutation c.3866-l-lG>A is a canonical splice site change and is 
predicted to remove the splice donor site of exon 5. An alternative 
spKce site is predicted at position -H99 in intron 5, if this alternative 
site is used it would lead to a frameshift and a premature stop 
codon (encoded by nucleotide 2-4 of intron 5). Both splice site 
mutations are thus predicted to lead to a frameshift and NMD 
could occur for at least part of the mRNA transcripts. The two 
other novel mutations, c.8222T>C (p.Phe2741Ser) and 
C.87670T (p.Arg2923*), found in two different families, are 
located in the region of the gene that encodes the tail region of 
MTOISA and are likely to cause a loss of function of this region. 
The p.Arg2923* mutation is present in the SH3 domain, which is 
involved in the protein -protein interactions [37]. Most of the 
previously reported mutations of MTOISA causing congenital 
severe to profound deafness were found in the motor head and the 
tail domains [6,38]. Collectively these results indicate that the 
motor head and tail regions of MTOISA are essential in the 
hearing process and any mutation in these regions is thus critical 
[39]. 

TMCl encoding a transmembrane protein is expressed in the 
neurosensory hair cells of the mouse cochlea [40] . The recurrent 
nonsense mutation in exon 7 of TMCl, c.lOOOT (p.Arg34*) is 
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the most common TMCl mutation reported in Pakistan. It leads to 
a truncated protein and leads to congenital severe to profound 
hearing loss [40,41,42,43]. Ben Said et al. [42] have previously 
described that the p.Arg34* is an old founder mutation found in 
several populations including the Pakistani population. It is likely 
that the family described here carries the mutation on the founder 
haplotype. Also a novel intronic mutation in the TMCl gene, 
c.362+18A>G was identified in the current study (Table 2), which 
creates a novel splice donor site and the insertion of 1 7 nucleotides 
as demonstrated in the minigene approach. 

BSND encodes the barttin protein, which is a vital subunit of the 
chloride and voltage-sensitive Ka (CLCNKA) and Kb (CLCNKB) 
channels in the inner ear and the kidney. Chloride channels and 
the barttin protein form heteromers which function in the 
recycling of the K""" ions in the inner ear and salt reabsorption in 
the kidneys [44,45,46,47,48,49]. As predicted by in silico analysis, 
the currendy identified BSND mutation p.Val33Leu residing in the 
transmembranal region might result in the malfunctioning of 
chloride channels CLCNKA and CLCNKB. However, in the 
current family there were no complaints of any renal problems. 

Deafness caused by mutations in the TMPRSS3 gene is bilateral 
and severe to profound with no defects of the middle ear and the 
vestibular system [50]. Mutations in TMPRSS3 can also cause 
progressive hearing loss with a postlingual onset [51]. TMPRSS3 
mutation C.7260G (p.Cys242Trp) in exon 8 affects the serine 
protease domain as predicted by 3D modeling. Previously another 
mutation in the same exon of TA4PRSS3, c. 6470T (p.Ar- 
g216Leu) has been shown to result in a failure of the protein to 
undergo proteolytic cleavage resulting in the inactivation of the 
sodium channel [52]. 

The homozygous variant p.Leu7Arg (c.20T>G) in MSRB3 
identified in the current study is located in the mitochondrial 
signal sequence and may result in mislocalization of the protein. 
Ahmed et al. [53] have shown the importance of the mitochon- 
drial isoforms when they found a mutation c.55T>C (p.Argl9*) to 
underlie hearing impairment. This mutation also resides in the 
signal sequence for mitochondrial localization. However, the 
variant c.20T>G (p.Leu7Arg) was also found in heterozygous 
state in 3 out of 89 (1.7%) ethnically matched controls. The 
mutation might be a founder mutation in the corresponding 
population. However, currently it is uncertain whether this variant 
is the cause of hearing impairment in the family and further studies 
are necessary for a definite conclusion on the pathogenic effect of 
the variant. 
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Conclusions 

In the present study 53% (16/30) of the families were found to 
carry causative mutations in GJB2 illustrating that the most 
frequently involved gene in deafness in the Pakistani population is 
GJB2 followed by MTOISA (13%, 4/30) and TMCl (10%, 3/30). 
Based on these results it is therefore suggested that as an initial step 
for the genetic diagnosis of deafness, GJB2 should be analyzed in 
Pakistani patients and if this gene is excluded then microsateUite 
markers flanking MY015A and TMCl should be genotyped to 
exclude those genes or to indicate mutation analysis. Homozy- 
gosity mapping is an effective approach to determine the mutated 
genes in consanguineous families. Although a large number of 
deafness genes and mutations have aheady been identified, our 
study demonstrates that the fuU mutation spectrum in these genes 
is stiU not defined. Identification of novel mutations in these genes 
is important for genetic counseling and can provide handles for 
further studies on protein function. Genetic counseling of the 
families is important to better inform couples about the risk of their 
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offspring to be hearing impaired. In addition, as can be seen in the 
current study, recurrent as well as novel mutations in known genes 
define the disease in the families which seem to have some 
population specificity. This information can be important in 

developing screening strategies. 

Supporting Information 

Figure SI Multiple-alignment of the corresponding 
stretches of protein sequences across different species. 

The blue color shading represents the intensity of conservation, 
where the dark blue shading represents highly conserved stretch 
while the light blue shading denotes moderate conservation of 
residues across different species. A) Amino acid sequence 
conservation of p.Gly200 across 9 species. B) Amino acid sequence 
conservation of p.Phe2741 across 11 species. C) Amino acid 
sequence conservation of p.Val33 across 7 species. D) Amino acid 
sequence conservation of p.Cys24 across 10 species. E) Non- 
conserved residue of p.Leu7 across 10 species. 
(TIF) 
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